Charge order, orbital order, and electron 
localization in the Magneli phase Ti407 



V. Eyert*, U. Schwingenschlogl, U. Eckern, 

Institut fiir Physik, Universitdt Augsburg, 86135 Augsburg, Germany 
Abstract 

The metal-insulator transition of the Magneli phase Ti407 is studied by means 
of augmented spherical wave (ASW) electronic structure calculations as based on 
density functional theory and the local density approximation. The results show 
that the metal-insulator transition arises from a complex interplay of charge order, 
orbital order, and singlet formation of those Ti 2>d states which mediate metal- 
metal bonding inside the four- atom chains characteristic of the material. Ti407 thus 
combines important aspects of Fe304 and VO2. While the charge ordering closely 
resembles that observed at the Verwey transition, the orbital order and singlet 
formation appear to be identical to the mechanisms driving the metal-insulator 
transition of vanadium dioxide. 
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1 Introduction 



Since tlie seminal work of Morin [1] the metal-insulator transitions of the 
early transition-metal oxides are motivating ongoing research. In particular, 
showing first-order transitions with a conductivity change of several orders 
of magnitude, VO2 and V2O3 are regarded as prototypical materials of this 
class. Recent theoretical work has revealed a strong influence via electron- 
phonon coupling of the structural degrees of freedom on the electronic prop- 
erties of VO2 and neighbouring rutile-type dioxides [2,3,4,5]. In contrast, the 
corundum- type sesquioxide has turned out to be subject to strong electronic 
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correlations [6,7], as seems to be the case also for Ti203 [8]. The relative im- 
portance of these mechanisms as well as their dependence on crystal structure 
and band filling are still unclear. 

In this situation the Magneli phases offer the opportunity of a more com- 
prehensive understanding [9]. These phases form a homologous series of the 
kind Mn02n-i (M = Ti, V) and are particularly suited for studying the differ- 
ences in crystal structures and electronic properties between the end members 
MO2 (n 00) and M2O3 (n = 2). This is due to the variation in d band 
occupation across the series as well as the fact that, as the general formula 
Mn02n-i = M2O3 + (n — 2)M02 suggests, the crystal structures can be viewed 
as rutile-type slabs of infinite extension and different thickness, separated by 
shear planes with a corundum-like atomic arrangement [9] . While in the rutile- 
type regions the characteristic metal-oxygen octahedra are coupled via edges, 
the shear planes have face-sharing octahedra as in the corundum structure. 
We have recently started a systematic investigation of the vanadium Magneli 
phases [10,11,12,13]; here we report on the first results for the titanium series. 

The Magneli phase Ti407 experiences a metal-insulator transition at 154 
K, which is followed by an insulator-insulator transition at 130 K with a 
thermal hysteresis of about 12 K [14,15,16]. Both transitions are first-order 
[15,16,17] and connected with a conductivity change of three orders of magni- 
tude each [14] . While a sharp decrease of the magnetic susceptibility has been 
reported for the 154 K transition, which has been interpreted as a change from 
Pauli paramagnetism to van Vleck-like behaviour [16], there is no consider- 
able change at lower temperatures [14,18]. In contrast, whereas the transition 
at 154 K goes along with increase in volume [17], a distortion of the crys- 
tal structure had been first observed only for the second transition [14]. From 
single-crystal x-ray diffraction data taken at 298, 140, and 120 K, Marezio et al. 
reported differentiation of the Ti-0 distances on going to the low-temperature 
phase, suggesting the formation of Ti^+ and Ti^"*" sites, which each occupy half 
of the characteristic four-atom chains and which contrast the average Ti'^'^+ 
valence observed at high temperatures [17,19,20]. In addition, the Ti"^"*" sites 
form short metal-metal bonds indicating singlet formation and explaining the 
changes in resistivity and magnetic susceptibility. At the same time the un- 
paired Ti^+ and, to a lesser degree, the Ti^+ atoms are displaced perpendicular 
to the chain axis off the center of the surrounding octahedron towards one of 
the oxygen atoms [17,19,20]. This pattern is similar to the situation found 
for Cr-doped VO2, where half of the chains pair and atoms on the remaining 
chains display a zigzag- like off-center displacement [3]. Finally, Marezio et al. 
reported large thermal displacements for the intermediate phase, suggesting 
that the structural changes are associated with the metal-insulator transition 
at 154 K but still lack long-range order [17,20]. The latter sets in at the lower 
transition, which is thus regarded as a disorder-order transition [17,20] as has 
been confirmed by specific-heat and EPR data [15,16]. These results led to 
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interprete the intermediate phase as a bipolaron hquid [15,16,21]. Yet, refined 
determination of the crystal structure revealed a fivefold superstructure for 
the intermediate phase, which shows the same structural characteristics as 
the low-temperature phase, and thus questioned understanding the conduc- 
tivity in terms of independent mobile bipolarons [22]. Recent photoemission 
and x-ray absorption measurements showed only minor effects of the insulator- 
insulator transition on the electronic properties but revealed drastic changes 
at 154 K [23,24]. The latter include significant down- and upshift of the oc- 
cupied and unoccupied Ti 3d states near the Fermi energy. In contrast, 
O 2p states extending from —9 to — 4eV remain essentially unaffected by the 
transition. 

Understanding of the physical properties of the Magneli phases is greatly fa- 
cilitated by an alternative unified representation of their crystal structures 
[10,11,12,13,17,20]. It is based on the notion of a hexagonal closed-packed ar- 
ray of oxygen atoms forming a regular 3D network of octahedra. Apart from 
a slightly different buckling, this network is the same for all Magneli phases 
including the rutile-type dioxide and the corundum-type sesquioxide. Differ- 
ences between the compounds Mn02n-i arise from filling the octahedra with 
metal atoms. Filled octahedra form chains of length n parallel to the pseu- 
dorutile Cprut axis, followed by n — 1 empty sites. While in the rutile structure 
these chains have infinite length, they comprise just two filled octahedra in 
the corundum structure, where they are perpendicular to the hexagonal c 
axis. Within the crystal, the metal chains are arranged in two different types 
of layers, which are interlaced by oxygen layers and alternate along ttprut- This 
gives rise to the aforementioned two types of chains, labelled (a) and (b). In 
Ti407, they comprise atoms Til/Ti3 and Ti2/Ti4, respectively (see Fig. 1 of 
Ref. [11], where this situation has been sketched for the isostructural V4O7). 
Due to alternation of the layers and relative shifts of the four-atom chains 
within the layers the end atoms Til and Ti2 are found on top of each other 
(note that in the notation by Marezio and coworkers [17,20] the chain center 
and end atoms are reversed). In the sesquioxide, these two atoms are usually 
designated as the c-axis pair. While neighbouring octahedra share faces along 
Oprut and feprut; mctal-metal bonding along all other directions within the lay- 
ers is via edges [3]. As a consequence, the atomic arrangement near the chain 
ends is corundum-like, whereas the chain centers correspond to the rutile-type 
regions. By virtue of the just sketched representation of the crystal structures 
it is possible to refer the symmetry components of the Ti 3d orbitals of all 
compounds to a common local coordinate system [3]. In this system the z- 
and of the local coordinate system are parallel to the apical axis of the 

local octahedron and the pseudorutile Cprut axis, respectively. 

In this letter we report on electronic structure calculations for the Magneli 
phase Ti407 using the crystal structures of both the room-temperature and the 
low-temperature phase. Our calculations reveal i) rather isotropic occupation 
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of the Ti 3d t2g states for the room-temperature structure, ii) significant elec- 
tron transfer from one half of the chains to the other in the low-temperature 
structure, iii) orbital order at the chains coming with strong metal-metal 
dimerization, and iv) strong electron localization, which leads the way to the 
observed metal-insulator transition. 



2 Methodology 



The calculations were performed using the scalar-relativistic augmented spher- 
ical wave (ASW) method [25,26]. In order to represent the correct shape of the 
crystal potential in the large voids of the open crystal structure, additional 
augmentation spheres were inserted. Optimal augmentation sphere positions 
as well as radii of all spheres were automatically generated by the sphere ge- 
ometry optimization (SGO) algorithm [27]. Self-consistency was achieved by 
an efficient algorithm for convergence acceleration [28]. Brillouin zone sam- 
pling was done using an increased number of k-points ranging from 108 to 
2048 points within the irreducible wedge. 



3 Results and Discussion 



Calculated partial densities of states (DOS) are displayed in Fig. 1. Here and 




Fig. 1. Partial DOS as resulting from the room-temperature (left) and 
low-temperature (right) crystal structure. Here and in the following figures slight 
broadening is due to the DOS calculation scheme [29]. 

in all following figures results for the room-temperature and low-temperature 
structure are given on the left and right, respectively. In Fig. 1, contributions 
from both the Ti 3(i and O 2p states are included. In addition, Ti contributions 
from the (a) and (b) chains containing Til/Ti3 and Ti2/Ti4, respectively. 
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are distinguished. All other orbitals play only a negligible role in the energy 
interval shown. 

Two groups of bands are identified, which extend from —0.5 to 2.0 eV and from 
2.4 to 4.6 eV and comprise 24 and 16 bands, respectively. They derive mainly 
from the Ti 3d states. In addition, 42 bands are found in the energy range 
from —8.5 to — 3.8eV, which trace back mainly to the O 2p states. However, 
p-d hybridization causes d and p contributions, respectively, below and above 
— 2eV reaching about 15% especially in the two upper groups. Crystal field 
splitting expected from the fact that the Ti atoms are located at the centers 
of slightly distorted octahedra is observed in the Ti 3d partial DOS. As a 
detailed analysis of the wave functions reveals, the lower group of Ti 3d derived 
bands arises almost exclusively from states of symmetry. In contrast, bands 
between 2.4 to 4.6 eV trace back to the Cg states. The energetical separation of 
the centers of gravity of these two groups of about 2.7 eV compares well with 
the value of 2.4 eV taken from the XAS measurements by Abbate et al. [23]. 
At the same time, the calculated position and width of the O 2p dominated 
group of bands is in very good agreement with the photoemision data of these 
authors. 

In general, these findings are very similar for the room-temperature and the 
low-temperature structure. In particular, the O 2p bands keep their energetical 
position and band width, again in agreement with the photoemssion data [23]. 
Differences show up on closer inspection of the Ti 3d dominated states. Most 
striking are the changes occurring in the occupied part of the Ti 3d t2g bands. 
While for the room-temperature structure the Ti 3d partial DOS in this en- 
ergy interval are almost the same for both chains, we observe strong increase 
and decrease of the Ti 3d partial DOS arising from the atoms of chain (a) and 
(b), respectively, for the low-temperature structure indicative of a consider- 
able charge transfer. The integrated difference amounts to ^ 0.8 electrons per 
formula unit. At the same time, states in the Ti 3d Cg group of bands, i.e. be- 
tween 2.4 and 4.6 eV, experience some rearrangement. In particular, on going 
from the room-temperature results to those obtained for the low-temperature 
structure, we find energetical down- and upshift, respectively, of the centers 
of gravity of these bands for chain (a) and (b). 

Both these changes can be attributed to the difference of the average Ti- 
O distances emerging in the low-temperature phase [17]. They lead to the 
assignment of Ti d^ and d^ charges contrasting the formal d^'^ valence of the 
room-temperature phase [17]. While in the latter phase the average Ti-0 bond 
length is the same for all titanium sites, in the low-temperature phase the av- 
erage bond length increases and decreases for chain (a) and (b) atoms, respec- 
tively. As a consequence, d electron charge is transferred from atoms Ti2/Ti4 
to Til/Ti3, leaving the d states in the (b) chains unoccupied. In addition, the 
shrinking and inflating of the oxygen octahedra strongly affect the bonding- 
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antibonding splitting especially of the cr-bonding states of Cg symmetry. This 
splitting decreases/increases in chains (a)/(b), leading to the observed down- 
/upshift of the Ti 3d dominated antibonding states. 

In order to investigate the charge redistribution coming with the transitions 
in more detail we display the Ti 3d partial DOS of all four titanium atoms 
separately in Figs. 2 and 3. In doing so we distinguish the three different 





Fig. 2. Partial 3d t2g DOS of (a) chain atoms Til and Ti3. 

d states contributing to this group of bands using the aforementioned local 
coordinate system as defined in Ref. [3]. In this system, the dx2_y2 states, 
pointing parallel to the pseudorutile Cprut axis, mediate cr-type metal-metal 
overlap within the chains. In contrast, the dyz orbitals take part in albeit 
weaker cx-type d-d overlap across the chains. Weak interchain coupling is also 
due to overlap of d^^^-y^ and dxz orbitals of neighbouring chains. 

For the high-temperature structure strong splitting of the d^2_y2 states of 
both Ti3 and Ti4 is observed. This is a consequence of the short Ti3-Ti3 and 
Ti4-Ti4 distances in the chain centers causing strong bonding-antibonding 
splitting of the dx2-y2 orbitals of these atoms. In contrast, distances between 
the chain center and end atoms are considerably larger, this leading to the 
somewhat unstructured shape of the dx2_y2 partial DOS of atoms Til and 
Ti2. Similar shapes are displayed by the dxy partial DOS of all atoms due 



6 



-1 -0.5 0.5 1 1.5 2 2.5 

(E - Ep) (eV) 



-1 -0.5 0.5 1 1.5 2 2.5 
(E - Ep) (eV) 




Fig. 3. Partial 3d t2g DOS of (b) chain atoms Ti2 and Ti4. 

to the very small metal-metal bonding these states take part in. Finally, the 
dyz states of all atoms show distinct splitting into bonding and antibonding 
states, reflecting overlap of these orbitals of atoms Ti3 and Ti4 (peaks at 0.2 
and 1.35 eV) as well as atoms Til and Ti2. In the latter case overlap is parallel 
to the pseudorutile a^rnt axis perpendicular to the chains across octahedral 
faces. This bond corresponds to the c-axis pairs in the sesquioxide. Indeed, 
the shape of the dyz partial DOS with its pronounced antibonding peak at 
1.6 eV and the tail of bonding states extending from —0.3 to 1.3 eV has been 
also observed for V2O3. 

On going to the low-temperature phase several striking changes occur. First, 
we observe the abovementioned increase and decrease of the partial DOS of 
the (a) and (b) chain atoms, respectively, in the occupied part of the spec- 
trum, which is a consequence of the charge ordering coming with the increased 
and reduced average Ti-0 bond-length as reported by Marezio and cowork- 
ers [17,20]. In the (b) chains smaller distances induce a stronger bonding- 
antibonding splitting of the hybridized Ti "id and O 2p states, hence, upshift 
of the antibonding d dominated bands. This effect is supported by the displace- 
ment of these atoms perpendicular to the chains. Eventually, we witness almost 
complete depopulation of the (b) chain states, i.e. Ti d^ valence. In contrast, 
on the (a) chains longer Ti-0 distances lead to reduced bonding-antibonding 
splitting, hence, energetical lowering of the d bands and the occupation reaches 
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nearly d^. 

Second, in addition to the charge ordering we find strong changes in the oc- 
cupation of the single d states of the (a) chain atoms. Due to the albeit small 
displacement of atoms Til and Ti3 perpendicular to the chain axis the d^z and 
dyz states of these atoms are subject to larger d-p hybridization and thus shift 
to higher energies. In contrast, both the Til and Ti3 d^^-y^ partial DOS show 
strong splitting into a double maxima structure with peaks near —0.25 and 
1.5 eV. This effect results from the dimerization of these two atoms observed 
below the metal-insulator transition, which leads to metal-metal overlap of 
the orbitals extending along the (a) chains and splitting of the corresponding 
bands into bonding and antibonding states. At the same, chain (b) atoms do 
not display any splitting of these states due to the uniform distance of more 
than 3.0 A between all atoms of this chain. 

The corresponding changes of the band structure as shown in Fig. 4 have dras- 
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Fig. 4. Electronic band structure. 

tic effects especially near the Fermi energy. While in the room-temperature 
structure bands in the lower part of the t2g group show strong dispersion lead- 
ing to the metallic behaviour, the low-temperature structure is characterized 
by two split-off bands. As expected from the partial DOS and as revealed by 
a detailed analysis of the wave functions, these weakly dispersing bands are 
almost exclusively of Til and Ti3 d^'^-y'^ character. Additional bands of the 
same symmetry, which likewise show a reduced dispersion, are found at the 
upper edge of the t2g group, supporting the previous interpretation in terms 
of bonding-antibonding splitting of these states. Yet, despite the separation of 
the lower d^2_y2 bands from all higher lying states there is still a semimetal- 
like overlap of about 0.1 eV, inhibiting the complete opening of the insulating 
band gap. 

In general, the situation resembles very much that found for the rutile-type 
compounds Ti02 and VO2, where the metal atoms are found in a d^ and d^ 
configuration, respectively. While the insulating titanate preserves the rutile 
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structure down to lowest temperatures maintaining especially the equidistant 
spacing of the metal atoms, vanadium dioxide shows distinct deviations from 
this structure at the metal-insulator transition. In particular, vanadium atoms 
experience pairing parallel to the 

'-'rut cl-^lS cLS well as zigzag-like antiferroelec- 
tric displacements perpendicular to this axis. In the calculation for VO2, these 
changes lead to splitting of the dx2-y2 bands as well as energetical upshift of 
the dxz and dyz states [3]. We may thus regard low-temperature Ti407 as com- 
posed of Ti02- and V02-like chains, which display the same behaviour as the 
respective dioxide. Common to VO2 and Ti407 is also the finding of a small 
semimetal-like overlap for the low-temperature phase, which we attribute to 
the shortcomings of the local density approximation. Yet, this weakness does 
not undermine the interpretation of the phase transitions of both compounds 
as arising from orbital order and strong bonding- ant ibonding splitting of the 
dx2_y2 bands due to the structural changes and strong electron-phonon cou- 
pling. In addition, Ti407 is characterized by distinct charge ordering, leading 
to separated d^ and d^ chains and laying ground for the efficacy of the afore- 
mentioned mechanisms. 



4 Conclusion 

In summary, according to electronic structure calculations for the Magneli 
phase Ti407, the metal-insulator transition of this material arises from the 
complex interplay of different types of ordering phenomena. While in the 
room-temperature structure all Ti "id t2g states display uniform occupation, 
the low-temperature structure is characterized by the formation of two types 
of titanium chains comprising Ti d^ and d^ states, respectively. This charge or- 
dering is accompanied by orbital ordering on the d^ chains and strong bonding- 
antibonding splitting of the respective d^^-y^ states coming with Ti-Ti dimer- 
ization and zigzag-type antiferroelectric displacement of these atoms in com- 
plete accordance with the findings for the d^ dioxide VO2. Finally, like in the 
latter compound, singlet formation on the d^i-yi states paves the way for the 
observed insulating behaviour. 

To conclude, Ti407 shares fundamental aspects with other cornerstone mate- 
rials of condensed matter physics. The charge ordering occuring at the metal- 
insulator transition is closely related to that discussed for the Verwey tran- 
sition of Fe304. In contrast, the orbital order and metal-metal dimerization 
accompanying the phase transition of Ti407 arise from the same mechanism 
as that driving the transitions of VO2 and Nb02, namely instability of the d^ 
atom chains towards the strong metal-metal bonding along the characteristic 
chains as well as antiferroelectric displacement perpendicular. This effect may 
thus be regarded as a universal mechanism common to many early rutile-type 
transition-metal compounds. 
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